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ABSTRACT 

We quantize the self-dual massive theory by using the Batalin-Tyutin Hamiltonian 
method, which systematically embeds second class constraint system into first class 
one in the extended phase space by introducing the new fields. Through this analysis 
we obtain simultaneously the Stiickelberg scalar term related to the explicit gauge- 
breaking effect and the new type of Wess-Zumino action related to the Chern-Simons 
term. 
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1 Introduction 



The Dirac method has been widely used in the Hamiltonian formalism [1] to quan- 
tize the second class constraint system. However, since the resulting Dirac brackets 
are generally field-dependent and nonlocal, and have a serious ordering problem be- 
tween field operators, these are under unfavorable circumstances in finding canonically 
conjugate pairs. On the other hand, the quantizations of the first class constraint 
systems [2,3] have been well appreciated in a gauge invariant manner preserving Becci- 
Rouet-Stora-Tyutin (BRST) symmetry [4,5]. If the second class constraint system can 
be converted into first class one in an extended phase space, we do not need to de- 
fine Dirac brackets and then the remaining quantization program follows the method 
of Ref. [2-5]. This procedure has been extensively studied by Batalin, Fradkin, and 
Tyutin [6,7] in the canonical formalism, and applied to various models [8-10] obtaining 
the Wess-Zumino (WZ) action [11,12]. 

Recently, Banerjee [13] has applied the Batalin- Tyutin (BT) Hamiltonian method 
[7] to the second class constraint system of the abelian Chern-Simons (CS) field theory 
[14-16], which yields first class constraint algebra in an extended phase space by in- 
troducing the new fields. As a result, he has obtained the new type of an abelian WZ 
action, which cannot be obtained in the usual path-integral framework. Very recently, 
we have quantized the nonabelian case, which yields the weakly involutive first class 
system originating from the second class one, by generalizing this BT formalism [17]. 
As shown in these works, the nature of second class constraint algebra originates from 
the symplectic structure of CS term, not due to the local gauge symmetry breaking. 
Banerjee, and Ghosh [18] have also considered a massive Maxwell theory, which has the 
explicit gauge-breaking term, in the BT approach. As a result, the extra field in this 
approach has identified with the Stuckelberg scalar. There are some other interesting 
examples in this approach [19]. 

In the present paper, we shall apply the BT Hamiltonian method [7] to the self-dual 
massive theory [20] revealing both the Stuckelberg effect [21] and CS effect [13,17] by 
using the BT quantization [7,18]. In section 2, since the BT formalism has developed 
quite recently, we first recapitulate this formalism by explicitly analyzing the pure CS 
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theory. In section 3, we apply the results discussed in section 2 to the well-known 
self-dual massive theory including the CS term in three dimensions, which is gauge 
non-invariant. By identifying the new fields p and A with the Stiickelberg scalar and 
WZ scalar, respectively, we obtain simultaneously the Stiickelberg scalar term related 
to the explicit gauge-breaking mass term and the new type of WZ action, which also 
includes the Stiickelberg scalar in order to maintain the gauge invariance related to the 
CS term. 



Now, we first recapitulate the BT formalism by analyzing the pure abelian CS 
model 



Since this action is invariant up to the total divergence under the gauge transformation 
5B^ = d^A, this action has a different origin of the second class constraint from the 
well-known massive Maxwell theory [18], which is due to the explicit gauge symmetry 
breaking term in the action. The origin of the second class constraints is due to the 
symplectic structure of the CS model. 

Following the usual Dirac's standard procedure [1], we find that there are three 
primary constraints 



2 The BT Formalism - The Pure CS Model 




(i) 



^0 = 



7T « 0, 

m + -me ij B j « (i = l,2) 



V t EE 



(2) 



and one secondary constraint, 



u 3 ee —mtijd l B 3 ~ 0, 



(3) 



obtained by conserving f2o with the total Hamiltonian, 




(4) 
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where H c is the canonical Hamiltonian, 

H c = J d 2 x [me ij B d i B j ], (5) 

and we denote x = (t,x) and two-space vector x = (x 1 ,^ 2 ) and e±2 = e 12 = 1 and 
Lagrange multipliers u°, u\ No further constraints are generated via this iterative 
procedure. We find that all rest constraints except tt = rr are superficially second 
class constraints. However, in order to extract out the true second class constraints, it 
is essential to redefine 003 by using Qi and ^2 as follows 

= d^i-^-me^BK (6) 

Then, f2 , f2 3 form the first class algebra, while f2 2 form the second class algebra 
as follows 

\j(x,y) = {toi(x),toj(y)} 

= ( _° m y^x-y); i,j = 1,2. (7) 

In order to convert this system into first class one, the first objective is to transform 
fli into the first class by extending the phase space. Following the BT approach [7], 
we introduce new auxiliary fields $ l to convert the second class constraint fi; into the 
first class one in the extended phase space, and assume that the Poisson algebra of the 
new fields is given by 

{&(x),&(y)} = u i '(x,y), (8) 

where u/- 7 is an antisymmetric matrix. Then, the modified constraint in the extended 
phase space is given by 

00 

tufa, fl", &) = fii + E ^ {n) ; ^ (n) ~ m n , (9) 

n=l 

satisfying the boundary condition, Cl^n^, B^, = 0) = The first order correction 
term in the infinite series [7] is given by 

^\x) = JdPyXijix^&iy), (10) 



and the first class constraint algebra of Qj requires the condition as follows 

Aij(a:, y)+f d 2 w d 2 z X ik (x, w)u kl (w, z)X lj (z, y) = 0. (11) 

As was emphasized in Ref. [13,17], there is a natural arbitrariness in choosing u/- 7 and 
Xij from Eq.(8) and Eq.(10), which corresponds to the canonical transformation in the 
extended phase space [6,7]. We take the simple solutions as 

w ij (x,y) = (J "q 1 ) S 2 (x-y), 

Xij(x,y) = J-^S^x-y), (12) 

which are compatible with Eq.(ll) as it should be. Using Eqs.(9), (10) and (12), the 
new set of constraints is found to be 

Cli = TTi + lme ij B j + (i = 1,2), (13) 

which are strongly involutive, 

{n a ,n p } = (a, (3 = 0,1, 2, 3) (14) 

with VLq = Q and f2 3 = Jl 3 . As a result, we have all first class constraints in the 
extended phase space by applying the BT formalism systematically. Observe further 
that only contributes in the series (9) defining the first class constraint. All higher 
order terms given by Eq.(9) vanish as a consequence of the choice Eq.(12). 

Next, we derive the corresponding involutive Hamiltonian in the extended phase 
space. It is given by the infinite series [7], 

oo 

H = H c + J2H (n) ; H {n) ~ ($*)", (15) 

n=l 

satisfying the initial condition, H^^B^,® 1 = 0) = H c . The general solution [7] for 
the involution of H is given by 

# (n) = -- J d 2 xd 2 yd 2 z &(x)uj tj (x, y)X^ k (y, z)g£~ 1] \z) , (n > 1), (16) 

where the generating functions are given by 

cf = ufUj, 

= {^\HW}o + {d?\HW}o (n>l), (17) 
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where the symbol O in Eq.(17) represents that the Poisson brackets are calculated 
among the original variables, i.e., O = (7^, B^). Here, and X l i are the inverse 
matrices of oj l i and Xy respectively. Explicit calculations yield, 

Gf ] = -mtij&B , (18) 

which is substituted in Eq.(16) to obtain 

#M = J cPxi-y/mid&^B ]. (19) 

In this case, there are no further iterative Hamiltonians. Thus, the total corresponding 
canonical Hamiltonian is 

H = H C + H {1 \ (20) 
which is involutive with the first class constraints, 

{n t ,H} = 0, (i = 1,2,3), 

{n ,H} = n 3 . (21) 

This completes the operatorial conversion of the original second class system with 
Hamiltonian H c and constraints Qi into first class with the involutive Hamiltonian 
H and constraints fij. Note that this Hamiltonian naturally generates the first class 
Gauss' law constraint f2 3 from the time evolution of Cl Q . 

Let us identify the new variables $ l as a canonically conjugate pair (A, ti\) in the 
Hamiltonian formalism, 

S 1 - -U A , 



jm 

^2 



$ 2 -> y/mX (22) 

satisfying Eqs.(8) and (12). Then, the starting phase space partition function is given 
by the Faddeev formula [3,22] as follows 

r 3 
Z = / VB^Vn^VWnx JJ S(n a )6(T p )det \ {Q a ,T p } \ e tS , (23) 

J a,/3=0 

where 

S = Jd 3 X (t^B" + 7T\\ - H) , (24) 



with Hamiltonian density Ti. corresponding to Hamiltonian H, which is now expressed 
in terms of (\,7T\) instead of $\ The gauge fixing conditions Tj are chosen so that 
the determinant occurring in the functional measure is nonvanishing. Moreover, Tj 
may be assumed to be independent of the momenta so that these are considered as 
Faddeev-Popov type gauge conditions [23]. 

We now perform the momentum integrations to obtain the configuration space 
partition function. The no, 7Ti, and 7r 2 integrations are trivially performed by exploiting 
the delta function 5(Q ) = 6(ir ), <S(f2i) = 5(tti + f B 2 + tt a ), and 5(Q 2 ) = 
S(7i2 — y-B 1 + m\), respectively. After exponentiating the remaining delta function 
5(Qs) = 5(—meijd l B^ + diir\ + md 2 X) with Fourier variable £ and transforming B° — > 
B° + ^, we finally obtain the action as follows 

S = J d 3 x[-^m€^ p B^d u BP + m\F 02 ], (25) 

where F 02 = d B 2 — d 2 B Q , and the corresponding Liouville measure just comprises the 
configuration space variables as follows 

3 

[Dfj] = VB i "VXD^5{Fqi + A) H{5(TplB + Z,B\\])}det \ {Q Q ,r^} |, (26) 

(3=0 

where 5(Fqi + A) is expressed by / T>i\\e~ % -I dZx ( F ° 1 +^) 7r >>. This action is invariant up to 
the total divergence under the gauge transformations as SB^ = <9 M A and <5A = 0. It 
is easily checked for consistency that starting from the Lagrangian (25) with a factor 
S(F 01 + A) in the measure part, one can exactly reproduce the set of all first class 
constraints Cl a and the involutive Hamiltonian (20). Note that we will show in the 
next section that the above 5-function, which is remaining in the measure part for the 
case of the pure CS theory, will be disappeared for the case of the non-pure CS theories 
[13,17] like the self-dual massive theory [20]. 



3 The Self-Dual Massive Model 

We consider the Abelian self-dual massive model [20] 

Ssd = J d 3 x [\m 2 B»B p - l -me, vp B^ '£?'], (27) 
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by using the useful results discussed in the previous section. Note that this action 
has an explicit mass term, which breaks the gauge symmetry as the case of the Proca 
model, and also the CS term, which has a different origin of the second class constraint 
system. Consequently, this action represents a second class constraint system, which 
can be easily confirmed by the standard constraint analysis. There are three primary 
constraints, 



7T « 0, 



~0, (i = l,2), 
and one secondary constraint, 

u 3 = m 2 B° - me ij d i B j « 0, 
which is obtained by conserving fl with the total Hamiltonian, 

H T = H C + J d 2 x[u°n + 11%}, 
where H c is the canonical Hamiltonian, 



d x 



m 2 {{B 1 ) 2 - {By} + me i3 B°d l Bi 



(28) 



(29) 



(30) 



(31) 



and u° and u l are Lagrange multipliers. No further constraints are generated via this 
iterative procedure. We find that all constraints are fully second class constraints. It 
is, however, essential to redefine u> 3 by using Qi and Q 2 as follows 



tt 3 = uj 3 + d i Vl i 
1 
2 



= 9V,, 



me ij d i B j + m 2 B°, 



(32) 



although the redefined constraints are still completely second class in contrast to the 
case of the pure CS theory. Otherwise, one will have a complicated constraint algebra 
including the derivative terms, which is difficult to handle. Then, the constraint algebra 
is given by 



{n a (x),n,3(y)} 

( 
0m 
-m 



\ m 







-m 2 \ 



/ 



S 2 (x- y y,a,P = 0,l,2,3, 



(33) 
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which reveals the simple second class nature of the constraints Q a (x). 

In order to convert this system into the first class one, the first objective is to 
transform Q a into the first class by extending the phase space. Following the BT 
approach [7], we introduce the matrix (8) as follows 



a/3. 



x,y) 



( 



V -i 



1 \ 
0-10 

1 
J 



5 2 {x-y). 



(34) 



Then the other matrix X a/3 in Eq.(10) is obtained by solving Eq.(ll) with A a/3 given 
by Eq.(33), 

/ m \ 



X a p{x,y) 






V o 




Tn 









m 

m J 



5\x-y). 



(35) 



There is an arbitrariness in choosing u a/3 , which would naturally be manifested in 
Eq.(34). This just corresponds to the canonical transformations in the extended phase 
space. However, as has also been evidenced in other calculations [13,17], this choice of 
Eqs.(34) and (35) brings about remarkable algebraic simplifications. 

Using Eqs.(9), (10) and (35), the new set of constraints is found to be 



Q = 7r + m$°, 



^3 



m + -me ij B j + s/m®\ 



(i = l,2), 



d^i - -me ij d i B j + m 2 B° + m$ 3 , 
2 



(36) 



which are strongly involutive, 

{ti a ,np} = 0. (37) 

Recall the $ M are the new variables satisfying the algebra (8) with u af3 given by Eq.(34). 
The next step is to obtain the involutive Hamiltonian. The generating functional 
is obtained from Eq.(17). It is noteworthy that there are only two terms VL a 
and i n the expansion (36) due to the intuitive choice (34) and (35). Explicit 
calculations yield, 



= m 2 B° 



m€ ij d l B j , 
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G 



(o) 



-m 2 B i - meij&B , 



= m 2 d t B\ 



which are substituted in Eq.(16) to obtain 



(38) 



H w = J d 2 x [m&diB* + mVme ij ^ i B j + y/m&diti - $ 3 (mB° - e ij d i B j 



(39) 



This is inserted back in Eq.(17) to deduce as follows 



G? 1 = sfmd,¥ + Hi* 3 , 



G« 
G?» 



(40) 



which then yield if( 2 ) from Eq.(16), 
1 



#(2) 



2 2 -v/m 2 



• (41) 



Since =0 (n > 2), the final expression for the desired involutive Hamiltonian 
after the n = 2 finite truncations is given by 



h = h c + hw + hv>, 



(42) 



which, by construction, is involutive, 



{n a ,H} = o. 



(43) 



This completes the operatorial conversion of the original second class system with 
Hamiltonian H c and constraints Q a into first class with Hamiltonian H and constraints 
Vt a . 

Before performing the momentum integrations to obtain the partition function in 
the configuration space, it seems appropriate to comment on the involutive Hamilto- 
nian. If we directly use the above Hamiltonian, we will finally obtain the non-local 
action corresponding to this Hamiltonian due to the existence of ^$ 3 9i$ 1 -term in 
the action when we carry out the functional integration over $ 1 or $ 3 later. Fur- 
thermore, if we use the above Hamiltonian, we can not also naturally generate the 
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first class Gauss' law constraint SI3 from the time evolution of the primary constraint 
Qq, which is the first class. Therefore, in order to avoid these problems, we use the 
equivalent first class Hamiltonian without any loss of generality, which differs from the 
involutive Hamiltonian (42) by adding a term proportional to the first class constraint 
Cl 3 as follows 

H' = H+ —<5> 3 tt 3 . (44) 
m 

Then, this Hamiltonian n consistently generates the Gauss' law constraint such that 
{Q ,H} = fi 3 . Note that non-locality may also be avoided by changing the order of 
performing the momentum integrals. But, in this case, one can not directly reproduce 
the original theory by fixing the unitary gauge as well as the Gauss' law constraint. 
Furthermore, when we act this modified Hamiltonian on physical states, the difference 
with H is trivial because such states are annihilated by the first class constraints. 
Similarly, the equations of motion for observable (i.e. gauge invariant variables) will 
also be unaffected by this difference since f2 3 can be regarded as the generator of the 
gauge transformations. 

We now unravel the correspondence of the Hamiltonian approach including both 
the Stuckelberg effect and the CS effect. The first step is to identify the new variables 
$ M as canonically conjugate pairs in the Hamiltonian formalism, 



;>° 
$1 



$ — > mp, 
1 



7TA, 



"m 

$ 3 - -7T„ (45) 
m 

satisfying Eqs.(8), (34) and (35). The starting phase space partition function is then 
given by the Faddeev formula, 

r 3 
Z= VB^Vn^VXVnxVpVnp JJ 5(fi a )5(r /3 )rfet | {Q a ,T p } \ e lS \ (46) 



a,/3=0 



where 



S' = J d 3 x (n^ + tt a A + n p p - H') (47) 



11 



with the Hamiltonian density H! corresponding to H', which is now expressed in terms 
of (p, tt p , A, nx) instead of <E> M . As in the previous section, the gauge fixing conditions 
r a may be assumed to be independent of the momenta so that these are considered as 
Faddeev-Popov type gauge conditions. 

Next, we perform the momentum integrations to obtain the configuration space 
partition function. The tt q , 7Ti, and 7r 2 integrations are trivially performed by exploiting 
the delta functions 5(Q ) = 5(n + m 2 p), 5(Qi) = 5(tti + |5 2 + tt\), and 
5(£l 2 ) = ^(^2 — y-B 1 + mA), respectively. After exponentiating the remaining delta 
function <5(f2 3 ) = S(—me i jd t B j + diir\ + m<9 2 A + m 2 B° + n p ) with Fourier variable £ 
as £(^3) = j'V^e^f^'^ 3 and transforming B° — > B° + £, we obtain the action as 



follows 



S = Jd 3 x{^m 2 B^B fl -^me^ p B^B' } 



+ p[-m 2 (B + £)- m 2 d i B i - \n?di&p + m 2 d 1 X - md 2 n x ] 
+ tt p [p - - f] + A[-m£ 2 + m 2 ^ 1 - m<9 2J B° - ^m 2 A] 

+ *x[\-B 1 -mB? + &&-±* x ]-±m 2 Z 2 }, (48) 

and the corresponding measure is given by 

3 

[Vp] = VB^VWirxVpVnpVZHiSiTplB + B\ \, p])}det \ {fJ a ,I>} |, (49) 

0=0 

where B° — > transformation is naturally understood in the gauge fixing condition 

Note that the original theory is easily reproduced in one line, if we choose the 
unitary gauge 

r = P , r 1 = 7r A , r 2 = A, r 3 = 7r p , (50) 

and integrate over £. Then, one can easily realize that the new fields <3> M are nothing but 
the gauge degrees of freedom, which can be removed by utilizing the gauge symmetry. 

Now, we perform the Gaussian integration over 7i p . Then all £ terms in the action 
are canceled out, and integrating over 7Ta, the resultant action is finally obtained as 
follows 

S = Sstuckelberg + S\VZ] 
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Sstuckewerg = J d 3 x{-m 2 (B^ + p - 2pd^) - -me^B^B"}, 
= J d*x{ l -m 2 {B, + 8,pf - m%(pB») - ^me^B^B?}, 
S W Z = J d 3 x{^[\ + Foi + m(5 2 + <9 2 p)] 2 

+ m\[F 02 - m{B 1 + d lP ) - ^mX]}, (51) 

where F 01 = d B 1 — diB , F 02 = d B 2 — d 2 B , and the corresponding Liouville measure 
just comprises the configuration space variables as follows 

[Dp] =VB»VXVpVtfl{5(r(ilB + Z,B\X,p])}det | {Q a ,T(,} | . (52) 

13=0 

This action S is invariant up to the total divergence under the gauge transformations 
as SB^ = d^A, 5p = —A, and 5X = 0. Starting from the Lagrangian (51) with the 
boundary term, we can easily reproduce the same set of all first class constraints Q a , 
and the Hamiltonian such that 

H = H c + j d 2 x ['K X diB -mTT X + m 2 pd i B i + mXd 2 B + m 2 XB 1 

+ ^n 2 -mn x d 2 p-^m 2 d t pd i p + m 2 Xd 1 p+^m 2 X 2 + ^n 2 p }. (53) 

Then, if we add a term proportional to the constraint Q3, i.e., — ^n p Cl 3 , which is 
trivial when acting on the physical Hilbert space, to the above Hamiltonian (53), we 
can obtain the original involutive Hamiltonian (42), which is canonically equivalent to 
Eq.(53). Furthermore, this difference is also trivial in the construction of the func- 
tional integral because the constraint ^3 is strongly implemented by the delta function 
8(Cl 3 ) in Eq.(46). Therefore, we have shown that the constraints and Hamiltonian fol- 
lowing from the Lagrangian (51) are effectively equivalent to the original Hamiltonian 
embedding. As results, through BT quantization procedure, we have found that the 
Stiickelberg scalar p is naturally introduced in the mass term, and this p as well as the 
WZ scalar A is also included in the new type of WZ action. 

Note that the gauge invariance of Swz should be maintained because the second- 
class constraint structure related to CS term only comes from the symplectic structure. 
We also note that the Wess-Zumino action in Eq.(51) is gauge invariant in spite of the 
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lack of the manifest Lorentz invariance. On the other hand, since the unitary gauge 
(50) recovers the manifestly Lorentz invariant original action, the actual invariance is 
maintained from the fact that the final result for the partition function Z is independent 
of the gauge fixing conditions. The local gauge symmetry of the Wess-Zumino action 
naturally also survives in the configuration space. This means that the origin of Swz is 
irrelevant to the conventional gauge-variant Wess-Zumino like action [8,10,11], which 
cancels the local gauge anomaly of the second-class system. Interestingly the choice 
of p = and A=0 does not recover the original theory in the Faddeev-Popov type 
gauges. Finally, it seems appropriate to comment on the action (51). If we ignore the 
boundary term in the Lagrangian (51), we cannot directly obtain the involutive first 
class Hamiltonian as the case of the Proca theory explained in Ref. [18] because the 
boundary term plays the important role in this procedure. 

In summary, we have recapitulated the Batalin-Tyutin method, which converts sec- 
ond class system into first class one, by analyzing the pure CS theory, which has the 
different origin of the second class structure. Then, we have applied this method to 
the Abelian self-dual massive theory including the CS term. As results, we have shown 
that if we ignore the boundary term in action (51), the direct connection with the usual 
Lagrangian embedding of Stiickelberg can be made by explicitly evaluating the momen- 
tum integrals in the extended phase space partition function using Faddeev-Popov-like 
gauges and identifying an extra field p introduced in our Hamiltonian formalism with 
the conventional Stiickelberg scalar. On the other hand, we have also obtained a 
new type of Wess-Zumino action containing the WZ scalar A, which also includes the 
Stiickelberg scalar p in order to maintain the gauge invariance of the Swz related to 
the CS effect in the action (51). 
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